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Abstract. Observations of low mass protostars which probe small enough size scales to be within likely CO 
depletion regions show the highest [NHhD]/[NH3] ratios yet measured, of 4-33%. These molecular D/H ratios are 
higher than those measured on larger scales, showing that deuterium fractionation increases towards protostellar 
cores. As in cold clouds, such high ratios can be produced by gas-phase ion-molecule chemistry in the presence of 
depletion. Grain surface chemistry is less likely to explain the deuterium enhancement, as it would require higher 
fractionation in ices than current models predict. The link between accretion, depletion and high molecular 
deuterium fractionation is strongly supported. 
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1. Introduction 

There are two main ways of producing molecular D/H 
enhancements in molecular clouds over the ~ 10 -5 ra- 
tio of [HD]/[H2]. Firstly, grain surface chemistry may en- 
hance molecular D/H ratios (Brown & Millar ll989all989bl 
Charnley et al. 11997(1 . Secondly, some key gas phase re- 
actions involving destruction of deuterated species run 
slower at low temperatures than the equivalent reactions 
with hydrogen, and this leads to molecular D/H enhance- 
ments where a cold gas phase chemistry has been active. 
Furthermore, in colder gas, depletion of heavy molecules 
such as CO results in an increase of [H 2 D + ]/[H;j~] and 
molecular D/H ratios (Brown & Millar 1989a, Roberts & 
Millar l2()()()al2T)00bl Rodgers & Charnley 11003) • 

Depletion of heavy molecules during freezeout in cold 
molecular clouds can result in molecular D /H ratios above 
10% (Roberts & Millar 2000), and can explain many of the 
high ratios observed in cold, dark clouds: DCO + fraction- 
ation in L1544 and L134N (Caselli et al. IH^OO! Tine et 
al. 12000)1 and ND 2 H in L134N (Roueff et al. 12000)) . 

Young protostars should show the effect of freezeout 
as well as cold dark clouds, though grain mantle release 
may also influence the ratios, as in IRAS 16293—2422 (van 
Dishoeck et al. I1995fl . A few studies have looked at molec- 
ular D/H ratios in low mass protostars: HDCO and DCN 
(Roberts et al. l2002)l : NH 2 D (Shah & Wootten l200Tl Saito 
et al. 12000)1 : DCO+ (Williams et al. H005)l . However, the 
70-90" beam size of most of these studies is a severe lim- 
itation as the largest D/H enhancements are expected on 
small scales, where densities are high. In L1544, heavy 
depletion occurs in a region 13000 AU across - only 90" 
diameter in nearby Taurus (Caselli et al. I1999f) . In order 
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to prove a general link between accretion and high D/H 
ratios, observations are needed of molecular D/H ratios 
in a sample of protostars with high enough resolution to 
probe the depletion regions. 

I present the results of such a survey of [NH 2 D]/[NH 3 ] 
ratios towards 11 protostars in Perseus, Taurus and Orion, 
observing NH 2 D with the IRAM 30m at 28" resolution 
(maximally 11000 AU at the distance of Orion) and cen- 
tering the observations on the protostellar dust peaks 
where the strongest depletion is likely to be found. 

2. Observations 

NH 2 D observations were made at the IRAM 30m on Pico 
Veleta in May 2000 and June 2002 towards 11 protostellar 
objects in Perseus, Taurus and Orion (Tabled . I observed 
the 85.93 GHz 1(1,1)-1(0,1) ortho transitions (6 hyperfme 
lines)(see. Olberg et al. 1985). The 1(1,1) state lies 20.7 K 
above ground. Integration times were 10 minutes (posi- 
tion switched) and 20 minutes for the two Orion sources; 
system temperatures were 165-190 K with typical sum- 
mer conditions of 4mm P WV (precipitable water vapour) . 
Observations of the two NGC 1333 sources in June 2002 
were made under 7mm PWV (T sys = 260-290 K). Noise 
levels were 50-90 mK on 40 kHz channels. Spectra were 
corrected for the beam efficiency B e g = 0.73. 

The main NH 2 D line was detected in all 11 sources. 
In the seven Perseus sources, the four outlying hyperfme 
transitions were also detected with peak brightness above 
5<7. In L1551IRS5, the hyperfme transitions were above 
3<T. In the other Taurus and Orion sources, only the main 
hyperfme line was detected. Spectra are shown in Fig. ^ 

NH3 observations at 24 GHz were made at the 
Effclsbcrg 100m in May 2002 using the 8192-channel 
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Table 1. Source positions and velocities 



Object 



RA (J2000) Dec (J2000) wlsr 
[ h m s 1 [ ' " 1 km s" 1 



Perseus 
















B5IRS1 


03 


47 


41.3 


+32 


51 


42 


10.2 


L1448mms 


03 


25 


38.8 


+30 


44 


05 


4.7 


L1448NW 


03 


25 


35.6 


+30 


45 


34 


4.7 


HH211 


03 


43 


56.8 


+32 


00 


50 


9.2 


IRAS03282+3035 


03 


31 


20.4 


+30 


45 


25 


7.1 


NGC1333 IRAS4A 


03 


26 


04.8 


+31 


03 


14 


7.0 


NGC1333 DCO+ 


03 


26 


06.6 


+31 


03 


08 


6.5 


Taurus 
















L1551IRS 


04 


31 


34.1 


+18 


08 


05 


6.2 


L1527 


04 


39 


53.5 


+26 


03 


05 


5.7 


Orion 
















RN043 


05 


32 


19.4 


+ 12 


49 


32 


10.4 


HH111 


05 


51 


46.3 


+02 


48 


30 


8.5 



spectrometer, which allows observation of the (J, K) = 
(1, l)-(4, 4) metastable transitions simultaneously. System 
temperatures were 40-50 K. Calibration was checked 
on W3(OH) and spectra were converted to Kelvin. The 
Effclsberg beamsize at 24 GHz is 37". NH 3 (1, 1) was de- 
tected in all sources; the NH 3 (1, 1) hyperfines and the (2,2) 
main line were detected in all sources except the two Orion 
sources (RN043 and HH111). The upper state energies for 
the (1,1) and (2,2) lines are 23.8 and 65.0 K respectively. 
The (3,3) and (4,4) lines were not detected. 



L1448'mms 




Fig. 1. NH2D (upper) and NH3 (lower) spectra. The tem- 
perature scale is Tmb in Kelvin. Note that the NH2D lines 
appear wider because of the factor 3.6 frequency difference 
between the NH3 and NH2D transitions. 



3. Analysis and results 

NH2D observations for the Perseus sources with well- 
detected hyperfines were reduced using CLASS method 
HFS to fit the six hyperfine lines (see Tine et al. (2000) for 
frequencies and line strengths) and derive optical depths 
and linewidths (Table |2J) • Where not all the hyperfines 
were detected, or the HFS fit gave too great an uncer- 
tainty in the optical depth, the main line was fitted with 
a single Gaussian (L1527, L1551IRS5, RN043, HH111). 

NH3(1,1) observations were reduced using CLASS 
method NH3(1,1) to fit the hyperfine structure and derive 
optical depths and linewidths (Table [5}. Excitation tem- 
peratures in the (1,1) line were 5-10 K assuming unity 
filling factor, but as the filling factor is uncertain these 
temperature estimates are of limited value. More use- 
fully, I derived the rotation temperature T12 using the 
method of Bachiller et al. (|1987j> . fitting the (2,2) main 
line with a single Gaussian. These temperatures are inde- 
pendent of filling factor and calibration. Rotation temper- 
atures were T12 = 11 — 14 K in the Taurus and Perseus 
sources. At these low temperatures, the rotational temper- 
ature is a good approximation to the kinetic temperature 
T12 ~ T kin (Walmsley & Ungerechts IT3£5f . For RN043 



in Orion, T12 = 14.6 K. The NH 3 (2,2) line in the other 
Orion source, HH111, was too weak for an estimate of T12. 

For the Perseus sources, optical depths were avail- 
able from hyperfine fits, and column densities and 
[NH2D]/[NH3] ratios were derived assuming LTE at 
T kin = Ti 2 . The resulting [NH 2 D]/[NH 3 ] ratios are 0.17- 
0.33 (Table El). 

For the Taurus sources L1527 and L1551 IRS5, op- 
tical depths were available for NH3 but not for NH2D. 
NH2D beam-averaged column densities were derived from 
the main hyperfine integrated intensities at Tki n = T12. 
The low [NH 2 D]/[NH 3 ] ratios, 0.04-0.05, likely reflect the 
fact that the NH2D column densities are beam-averaged 
whereas the NH3 optical depths refer to higher column 
density material on smaller scales, so the ratios are lower 
limits. 

For the Orion sources RN043 and HH111, NH3(1, 1) 
optical depths were much less than one and very uncer- 
tain, so the NH3 beam averaged column density was cal- 
culated from the integrated intensity of the main hyper- 
fine line. NH2D and NH3 column densities were calculated 
at 14.6 K (the measured T 12 in RN043). The ratios are 
additionally divided by the ratio of the beam areas to 
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Table 2. NH2D and NH3 observed line parameters. Uncertainties (in brackets) come from either formal errors on 
optical depth from the hyperfine fits or 30% absolute calibration uncertainties where column density was derived from 
integrated intensity. 



Sonrcp 

LJ V ) III V . 






NH 2 D 85.93 GHz 










NKUfl 1) 






NH^2 2) 








J I MB UjV 


Af 




7" 




J A MB uu 


Av 




J i MB uu 








K km s -1 


km s 


-1 






K km s _1 


km t 




K km s" 1 


R5TRS1 


3.22 


(0 74) 




0.46 


(0.03) 


3.74 


(0.06) 




0.58 


(0.00) 




T ,1 448mm« 


2.43 


(0.29) 




0.88 


(0.02) 


2.79 


(0.04) 




1.00 


(0.00) 


1 93 (0 05) 


L1448NW 


3.66 


(0 12) 




0.83 


fo 01) 


5.37 


fO 04) 




0.89 


(0.00) 


1.41 (0.03) 


HH211 


4.54 


(0.28) 




0.49 


(0.01) 


2.70 


(0.15) 




0.68 


(0.01) 


0.53 (0.03) 


IRAS03282 


3.39 


(0.18) 




0.50 


(0.01) 


3.53 


(0.15) 




0.59 


(0.01) 


0.39 (0.02) 


NGC1333 IRAS4A 


1.46 


(0.29) 




1.38 


(0.04) 


1.66 


(0.50) 




1.28 


(0.06) 


0.78 (0.15) 


NGC1333 DCO+ 


1.84 


(0.25) 




1.19 


(0.02) 


2.54 


(0.08) 




1.20 


(0.01) 


0.95 (0.04) 


L1527 






0.247 (0.074) 


0.83 


(0.13) 


1.51 


(0.01) 




0.50 


(0.01) 


0.15 (0.02) 


L1551IRS5 






0.680 (0.204) 


0.60 


(0.03) 


2.14 


(0.02) 




0.66 


(0.01) 


0.52 (0.03) 


RN043 






0.301 (0.090) 


0.68 


(0.07) 






0.791 (0.237) 


1.04 


(0.04) 


0.18 (0.02) 


HH111 






0.177 (0.053) 


0.51 


(0.07) 






0.964 (0.289) 


1.03 


(0.04) 





Table 3. 

[NH 2 D]/[NH 3 ; 



Temperatures, 
ratio. 



column densities and 



Source 


T12 


iV(NH 2 D) 


iV(NH 3 ) 


[NH 2 D]/ 




K 


xlO 14 cm" 2 


xlO 14 cm" 2 


[NH 3 ] 


R5IRS1 


10.8 


2.12(0.49) 


13.4(0.2) 


0.18(0.04) 


L1448mms 


13.7 


3.95(0.47) 


18.0(0.3) 


0.20(0.02) 


L1448NW 


11.7 


4.74(0.16) 


29.7(0.2) 


0.17(0.01) 


HH211 


12.8 


3.85(0.24) 


11.6(0.6) 


0.33(0.03) 


IRAS03282 


11.3 


2.52(0.13) 


12.9(0.6) 


0.22(0.01) 


NGC1333: 










IRAS4A 


14.2 


3.90(0.78) 


13.8(4.2) 


0.25(0.09) 


DCO+ 


12.8 


3.76(0.51) 


19.2(0.6) 


0.19(0.03) 


L1527 


10.8 


0.16(0.05) 


3.79(0.02) 


0.04(0.02) 


L1551IRS5 


12.6 


0.42(0.12) 


8.88(0.07) 


0.05(0.01) 


RN043 


14.6 


0.18(0.05) 


0.46(0.14) 


0.22(0.10) 


HH111 




0.10(0.03) 


0.56(0.17) 


0.11(0.06) 



give [NH 2 D]/[NH 3 ] ratios of 11-22%. This assumes that 
the source is centrally condensed and N(NH3) is higher 
on the NH2D beam scales. The correction is conservative, 
and the [NH 2 D]/[NH 3 ] ratios for the Orion sources could 
be a factor of 1.75 higher if the NH 2 D is extended, or 
19-39%. 

When calculated without optical depth information 
the ratios are less certain, as the absolute calibration un- 
certainty at each frequency is ~ 30%. In contrast, where 
optical depths are measured the [NH 2 D]/[NH 3 ] ratio does 
not depend on the absolute calibration. All the column 
density determinations are of course limited by the lack 
of information about source structure from single spec- 
tra. Foreground self-absorption and multiple components 
(which particularly affect NGC1333) and the beam size 
difference between Effelsberg and the 30m cannot be fully 
accounted for with the present data. 

4. Discussion 

In the Perseus sources where the ratios are best deter- 
mined, [NH 2 D]/[NH 3 ] = 17-33%. The highest ratios de- 



rived by Shah & Wootten ©001) and Saito et al. (2000) 
are 13% for sources also on this sourcelist. Differences in 
excitation assumptions may also have an effect, but the 
lower ratios derived by Shah & Wootten (j2001|l can eas- 
ily be explained if the fractionation is enhanced on small 
scales, towards the protostellar cores, as their beam is 
larger (80-90" compared to 27"). Saito et al. pToTTjl 's 18" 
resolution observation towards Bl is 15" away from the 
dust peak (Matthews & Wilson 2002) and shows a lower 
ratio, which also supports the idea that [NH 2 D]/[NH 3 ] re- 
duces away from the cores. (No conclusion can be drawn 
from their uncertain NH 2 D column in L1448). The ratios 
for protostars are also higher than those towards L134N, 
presumably because the protostars are denser and more 
depleted ([NH 2 D]/[NH 3 ]= 0.1; Tine et ab ETTO]) . 

Why should [NH 2 D]/[NH 3 ] be higher close to the pro- 
tostars? Is cold gas-phase chemistry in the presence of 
freezeout enhancing ratios in the cores, or are highly frac- 
tionated grain ices evaporating near the young stars? 

Chemical models indicate that [NH 2 D]/[NH 3 ] ratios 
of 30% can be produced by ion-molecule chemistry if 
heavy molecules are depleted (see Fig. 3 of Roberts & 
Millar EOOnS Fig. 2 of Rodgers & Charnley IWTTl . The 
ratios can be even higher if branching ratios for disso- 
ciative recombination favour deuterium retention (Lis et 
al. 2002). The corresponding fractionation required in 
molecular ions is similar or less than that already mea- 
sured in L134N (Rodgers & Charnlev I21KJT1 N 2 D+ obser- 
vations by Tine et al. I2000[) . Gas-phase chemistry with 
freezeout and depletion can therefore relatively easily ex- 
plain the high [NH 2 D]/[NH 3 ] ratios. 

An alternative hypothesis is that removal of highly 
fractionated grain mantles enhances ratios. Infrared obser- 
vations suggest that NH 3 may be present in ices, though 
the level is still uncertain (Lacy et al. [T998, Gibb et 
al. IMiUl Dartois & D'Hendecourt HSU), but NH 3 ice 
alone evaporates at 60-70 K; in a water matrix tempera- 
tures greater than 90 K are required. The temperatures of 
our sources, which are mostly Class protostars, rule out 
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much thermal ice evaporation. Only B5IRS1 and L1551 
have bolometric (dust) temperatures above 60 K. The 
source with the highest [NH 2 D]/[NH 3 ] ratio, HH211, is 
probably extremely young (Gueth & Guilloteau [1223 and 
has a bolometric temperature of only 30 K. Kinetic tem- 
perature estimates for the NH3 are all less than 15 K. 

The possibility remains that highly fractionated grain 
mantles are being removed by non-thermal evaporation, 
ie. slow, non-dissociative shocks. [NH2D]/[NH3] ratios 
are around 10% in regions where grain evaporation is 
known to be occuring (hot cores and IRAS 16293; van 
Dishoeck et al. 1995), and may originate in the ices, as may 
high [D 2 CO]/[H 2 CO] ratios (Ceccarelli et al. 1998,2002). 
Enhancements of another molecule, CH 3 OH, in a simi- 
lar sample of cores, are also attributed to non-thermal ice 
evaporation (Buckle & Fuller [2002). Two conditions must 
be satisfied if a non-thermal removal mechanism operates. 
Firstly, it must reproduce the low linewidths of the NH3 
(< 1 km s _1 ) and NH2D (linewidths the same or lower), 
which suggest that NH 2 D is tracing quiescent rather than 
shocked material. Secondly, if the release of ice material is 
to explain all the NH 2 D observed, the [NH 2 D]/[NH 3 ] ratio 
in the ices must be as high or higher than that measured 
in the gas phase, because of dilution by low [NH 2 D]/[NH3] 
ratios in regions of the core without ice evaporation. 

Observations suggest that, in these sources, ices can 
only be released in a small fraction of the core, otherwise 
CH3OH enhancements of several hundred percent would 
be expected, as seen in outflows and hot cores, rather than 
the moderate 10% observed (Buckle & Fuller f2T)(l2*jl . and 
column densities of NH3 would be enhanced by orders of 
magnitude over the surrounding molecular clouds, which 
is not the case (Bachiller et al. IM7I ITMT1 Benson 
& Myers Anglada et al. lDSUjl . 

To obtain the observed ratios even without dilution, 
the statistical treatment of Rodgers & Charnley ( 2001 ) fol- 
lowing Brown & Millar I|1989all989b|l requires an atomic 
D/H ratio in the gas phase at the time of freezeout of 
0.8-1.6. This is at the limit of what either steady-state 
or depletion models predict, as these cannot produce 
atomic D/H > 0.1 (Roberts & Millar 2000a; Roberts, priv. 
comm.). 

With dilution, therefore, current theories of grain man- 
tle chemistry cannot explain the NH 2 D enhancement, but 
would require either higher [D]/[H] ratios during freeze- 
out or preferential formation of deuterated molecules in 
the ice. 

The [NH 2 D]/[NH 3 ] ratios of - 0.2 in NGC1333 is 
consistent with ND3 column densities (van der Tak 
et al. I2HU21 Lis et al. l2*U02j l. using the models of 
Rodgers & Charnley i|2001[l . This is in contrast to earlier 
[NH 2 D]/[NH 3 ] ratios of 10% which predicted [ND 3 ]/[NH 3 ] 
an order of magnitude lower than measured. As the theory 
predicts [ND 3 ]/[NH 3 ]oc[NH 2 D]/[NH 3 ] 3 , ND 3 predictions 
are very sensitive to [NH 2 D]/[NH3] ratios, and more obser- 
vations of both isotopomers are needed to say if the theo- 
ries for multiple deuteration arc correct. Measurements of 
ND 2 H in the same sources would provide a further test: 



[ND 2 H]/[NH3] ~ 0.03 is expected from the models, for 
the given [NH 2 D]/[NH3] ratios. Detailed maps of individ- 
ual sources would also test the relationship between high 
deuteration and CO depletion. 

5. Summary and conclusion 

[NH 2 D]/[NH3] ratios in low mass protostellar cores on 
10,000 AU scales can reach 30% percent. Such high ra- 
tios can be simply explained by recent models of cold gas- 
phase chemistry with depletion due to freezeout (Roberts 
& Millar 2000a, 2000b|). Current models of ice formation 
even under conditions of high depletion have difficulty ex- 
plaining the high fractionation observed. 
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